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ABSTRACT 

We have discovered a binary AGN in the galaxy Mrk 739 using Chandra and Swift BAT. We 
find two luminous (L2-10 keV=l-l x 10 43 and 1.0 x 10 42 erg s _1 ), unresolved nuclei with a projected 
separation of 3.4 kpc (5.8±0.1") coincident with two bulge components in the optical image. The 
western X-ray source (Mrk 739W) is highly variable (x2.5) during the 4-hour Chandra observation 
and has a very hard spectrum consistent with an AGN. While the eastern component was already 
known to be an AGN based on the presence of broad optical recombination lines, Mrk 739W shows no 
evidence of being an AGN in optical, UV, and radio observations, suggesting the critical importance of 
high spatial resolution hard X-ray observations (>2 keV) in finding these binary AGN. A high level of 
star formation combined with a very low L[ O ni]/^2-i0 keV ratio cause the AGN to be missed in optical 
observations. 12 CO observations of the (3-2) and (2-1) lines indicate large amounts of molecular gas 
in the system that could be driven towards the black holes during the violent galaxy collision and be 
key to fueling the binary AGN. Mrk 739E has a high Eddington ratio of 0.71 and a small black hole 
(log Mbh=7.05±0.3) consistent with an efficiently accreting AGN. Other than NGC 6240, this stands 
as the nearest case of a binary AGN discovered to date. 

Subject headings: galaxies: active - galaxies: individual: Mrk 739 - galaxies: interactions - X-rays: 
galaxies 



1. INTRODUCTION 

The detection and frequency of binary AGN provide 
constraints on models of galaxy formation and an im- 
portant test of the merger-driven AGN model. Hierar- 
chical merger models of galaxy formation and the fact 
that almost all massive galaxies have supermassive black 
holes suggest tha t binary black h oles should be common 
in galaxies (jVolonteri et al.ir2Q03h . If galaxy mergers are 
the prime way to 'ignite' the central source by sending a 
large amo unt of gas into the center region and triggering 
the AGN (|di Matteo et alJl2QQ5h , then we expect some 
fraction of these binary black holes to be actively growing 
simultaneously, thus creating a binary AGN. 

Despite their theoretical importance, only a handful 
of close binary AGN (<5 kpc projected separation) have 
been discovered. The two clearest cases are the Chandra- 
detected double nucleus in the lu minous infrared galaxy 
(LIRG; Li R >lO n L ) NGC 6240 (jKomossa et alJl2003[ ) 
with a projected separation of 1 kpc at a distance of 
103 Mpc and the LIRG Mrk 463 at 3.8 kpc separation 
and a distance of 220 Mpc. Recently some likely bi- 
nary AGN have been discove red based on dou ble-peaked 
[O III] A5007 emission lines (jLiu et al.ll20l(l . However, 
there is still some question whether these systems are 
'true' binary AGN or are a single AGN with an asymmet- 
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ric distribution of outflowing gas in the n arrow line region 
([Smith et all 120101 : iFischer et all I2011D . Unfortunately, 
these systems are at higher redshifts with extremely close 
separations where the resolution of Chandra is unable to 
resolve these objects to confirm their binary AGN nature. 

As part of our Chandra program o f following up th e 
close mergers detected by Swift BAT (Ko ss et aLl l2010). 
we have discovered a binary AGN in Mrk 739 with a 3.4 
kpc separation at a distance of 130 Mpc. The binary 
AGN is particularly interesting because it shows no ev- 
idence of being an AGN in the optical, UV, or radio. 
Other than NGC 6240, this stands as the nearest case of 
a binary AGN discovered to date. 

2. OBSERVATIONS AND DATA ANALYSIS 

In the following subsections, we describe the observa- 
tions and analysis of Mrk 739. Throughout this work, we 
adopt tt m = 0.3, Q A = 0.7, and H = 70 km s -1 Mpc -1 
to determine distances. At the redshift of Mrk 739, 1" 
corresponds to 580 pc. 

2.1. Optical: SDSS Imaging and Gemini Optical 
Spectroscopy 

Mrk 739 was imaged by the SDSS on March 10, 2005. 
Using a Sersic profile with a fixed bulge (n=4), we fit 
the optical nuclei usin g two-dimensiona l surface bright- 
ness fitting fGALFIT: 1Peng et al.ll2002h . In Mrk 739E, 
a point source component was used to measure the AGN 
light since it has a broad line region (BLR). 

We observed Mrk 739 with Gemini on February 7, 
2011. Both nuclei were observed simultaneously in the 
B600-G5307 grating with a 1" slit in the 4300-7300 A 
wavelengt h range. The expo sure totaled 37 minutes. 
We follow IWinter et aTl (|2010h for correcting Milky Way 
reddening, starlight continuum subtraction, and fitting 
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AGN diagnostic lines. To correct our line ratios for ex- 
tinction, we use the narrow B aimer line ratio (Ha/H/3) 
assuming an intrinsic ratio of 3.1 and the iCardelli et all 
(1989) reddening curve. 

2.2. UV and Radio: XMM -Newton and VLA 

Mrk 739 was observed in the UV with XMM-Newton 
in June 2009 (PI Brandt, Vasudevan in prep.). We follow 
the XMM ABC guide for photometry. We also analyzed 
archival VLA observations with times of 33 minutes at 
1.49 GHz and 38 minutes at 4.86 GHz. 

2.3. X-rays: Chandra 

Chandra observed Mrk 739 on April 22, 2011 with an 
exposure time totaling 13 ks. Sub pixel event reposition- 
ing was applied to improve the resolution of the image. 
Two extraction regions of 1.5" radius were used for spec- 
tral fitting and timing analysis with CI AO version 4.3. 
To fit the X-ray spect ra, we used a fixed Galactic pho- 
toelectric absorption (jKalberla et al.l l2QQ5h . a floating 
photoelectric absorption component at z=0.0297, and a 
power law. For the eastern source (Mrk 739E), we also 
include a pileup model because mild pileup (10%-20%) 
is expected based on pixel count rates. 

2.4. Submillimeter: CO Observations 

The 12 CO (2-1) and (3-2) molecular lines of Mrk 739 
were observed with the James Clerk Maxwell Telescope 
(JCMT) on March 12-13, 2011. The A3 (211-279 GHz) 
and HARP (325-375 GHz) receivers were used. The spec- 
tra were co-added, binned, and fitted with linear base- 
lines. To calculate velocity-integrated line flux densities, 
we assumed an aperture efficiency of 0.60 and 0.53 for th e 
A3 and HARP receiver and followed iGreve et al.l ()2QQ9h . 

3. RESULTS 

3.1. Optical Properties of Mrk 739 

Two hard X-ray point sources coincide with the best 
fit model of the optical light from the bulge components 
(Figure 1). The bulge magnitudes are m r =14.03±0.15 
for Mrk 739E and ra r =13.75±0.15 for Mrk 739W. The 
small difference in apparent magnitudes suggests a major 
merger between the two galaxies. 

The optical spectra for both sources are shown in Fig- 
ure 2. Mrk 739E shows broad lines (FWHM H/3=2960 
km s _1 and Ha=2120 km s _1 ) consistent with a Seyfert 
1. Mrk 739E also has strong [Fe VII] 5721 and 6087 A 
emission, a feature of some Seyfert 1 galaxies indicative 
of highly ion ized material near the central AGN (e.g., 
lVeilleuxl [r988). In Mrk 739W, there are narrow lines at 
the spectral resolution of instrument (FWHM=280 km 
s _1 ). Based on the Balmer decrement, E(B — V)=0.26 
for Mrk 739E and E(B - V)=0A3 for Mrk 739W. Mrk 
739W is classified as a starburst using the [OI] /Ha and 
[SII]/Ha line diagnostics and a composite ga laxy using 
the [Nil] /Ha diagnostic dKewlev et al.l l2006f ). Assum- 
ing that all of the Ha lum inosity is from star formation 
and using iKennicuttl (jl998l ). the estimated star formation 
rate (SFR)=0.3 M yr" 1 . 

3.2. UV, Far- Infrared, and Radio 




Figure 1. Top: composite SDSS gri filter image of Mrk 
739. Bottom: SDSS r-band image overlaid with contours from 
Chandra (blue) at same scale. The X-ray sources are coincident 
with the bulge components seen in the optical for Mrk 739E and 
Mrk 739W. 

The UV image from XMM-Newton Optical Mon- 
itor shows sources coincident with the hard X-ray 
sources (Figure 3). In the UV, muvwi=15.4±0.1 and 
^UVM2=16.6±0.1 for Mrk 739E and muvwi=16.6±0.1 
and ra UV M2=16.4±0.1 for Mrk 739W. The spectral in- 
dex connecting 2500 A and 2 keV, a x, is - 1.10±0.04 for 
Mrk 739E, typical of an AGN (-1.15±0.24: ISteffen et all 
I2006f ). In Mrk 739W, a O x=-1.53±0.1, suggesting the UV 
is dominated by star formation. Assuming all the UV 
emis sion in Mr k 739W is from star formation and follow- 
ing [Kennicut3 ([ 19981 ). we estimate SFR=0.6 M yr _1 . 

Mrk 739 was detected in all bands by IRAS. The mea- 
sured luminosity is log L 0) tr=1O.9 and log L 0? fir=1O.6. 
Following the Kennicutt (1998) relationship between far- 
infrared (FIR) luminosity and star formation rate, we 
estimate SFR=6.9 M yr _1 . 

Mrk 739E was detected at both 1.49 and 4.86 GHz 
(Figure 3), with a flux density of 2.6±0.2 mJy at 1.49 
GHz and 0.5±0.2 at 4.86 GHz (after convolving to match 
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Figure 2. Upper panel: optical spectra of Mrk 739E and Mrk 739W. Lower panel: emission line ratios of Mrk 739W, the narrow line 
source discovered to be an AGN in X-rays. Mrk 739W is consistent with star formation in the [OIJ/Hai and [SIIJ/Ha diagnostics and a 
composite galaxy in the [Nil] /Ha diagnostic. 



the 1.49 GHz resolution). This is a spectral index (S u = 
K x v~ a ) of a =1.2±0.5. 

The VLA data for Mrk 739W are consistent with re- 
solved star formation and show no signs of an AGN. 
The emission at 1.49 GHz has an integrated flux den- 
sity of 2.6±0.2 mJy. The 4.86 GHz convolved data also 
shows resolved emission with an integrated flux density 
of 1.3±0.3 mJy. The spectral index for Mrk 739 W is a = 
0.8±0.3, consistent with optically thin synchrotron emis- 
sio n from supernova e found in star forming galaxies. We 
use lYun et al.l (|2QQll ) to convert the 1.49 GHz luminosity 
of Mrk 739W to a SFR of 3.1 M© yr" 1 . 

3.3. X-ray 

We detect two hard X-ray sources coincident with the 
eastern optical nucleus (Mrk 739E) and western nucleus 
(Mrk 739W). Both sources show hard X-ray spectra ex- 
tending out to 10 keV (Figure 4). In the 2-10 keV band, 
we find a FWHM of 0.48" ± 0.05 (280 pc) for Mrk 739E 



and 0.51" ± 0.07 (295 pc) for Mrk 739W. The Chandra 
spectra of Mrk 739E is well fit (x 2 =1.4) by a power 
law (F oc £T r+2 ) with photon index r = 2.1 ± 0.1 
and A^h = 1.5 ± 0.2 x 10 21 cm -2 consistent with a 
Seyfert 1. Mrk 739W is well fit (C-stat/dof=0.8) by 
a harder power law with more absorption and a pho- 
ton index r = 1.0 ± 0.2 and Ah = 4.6 ± 0.1 x 10 21 
cm -2 . While positive residuals do exist at the location 
of the neutral 6.4 keV iron Ka line in Mrk 739 W, there 
are too few counts to confirm its existence. The 2-10 
keV absorption-corrected luminosities are L2-iokeV = 
1.1 xlO 43 and 1.0 x 10 42 erg s -1 for Mrk 739E and Mrk 
739W, respectively. An archival XMM-Newton 2009 
observation of Mrk 739 is unable to resolve the emis- 
sion to either source, but shows L2-10 keV=l-0 x 10 43 and 
r = 1.92 ± 0.02 consistent with the Chandra spectra of 
Mrk 0739E. 

Timing analysis is a critical part of identifying AGN 
since nearly all show variability. Bins of 1.5 ks were cho- 
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Figure 3. Top: XMM-Newton UVM2 image of Mrk 739. Bot- 
tom: the grey-scale image is the 1.49 GHz VLA data, while the 
contours are the 4.86 GHz data convolved to 1.49 GHz beam. 
White x's indicate the Chandra hard X-ray positions. Both the 
UV and radio data in Mrk 739W show extended emission consistent 
with star formation. 

sen to ensure >20 counts per bin (Figure 4). We find 
statistically significant variations in the fluxes of both 
sources. There is a factor of 2.5 change in flux for Mrk 
739W and 0.3 for Mrk 739E during the ^4-hour obser- 
vation. 

Mrk 739 was detected by Swift BAT with 
£14-195 kev=2.4±0.5 x 10 43 erg s~\ T=2.6 ± 0.4, and 
a factor of two variability over five years (Figure 4). Be- 
cause of the variability and much steeper photon index 
of Mrk 739W, it is difficult to identify the source of the 
£14-195 kev luminosity (Mrk 739E, Mrk 739W, or both). 

3.4. Detection of a Binary AGN 

The hard X-ray (>2 keV) band provides one of the best 
tools for finding AGN since it is less affected by contami- 
nation and absorption and can only be produced in large 
amounts by AGN. Our discussion of this binary AGN will 
be limited to Mrk 739 W since Mrk 739E was already 
known to be an AGN based on its BLR (Ne tzer et al.l 



[19871) . 

In Mrk 739W, the hard X-ray emission is point-like 
at the location of one of the bulge components in the 
galaxy consistent with an AGN. Energetic phenomena 
related to vigorous star formation such as OB stars, X- 
ray binaries, and SN shocks produce hard X-rays. How- 
ever, in star forming regions the dominant X-ray emission 
is from point-like u ltraluminous X-ray sources (ULXs; 
iBertram et al.ll2007[ ). ULXs are by definition not located 
at the centers of galaxies where the central supermassive 
black hole resides. 

The luminous hard X-ray emission, hard photon in- 
dex, and time variability of Mrk 739W also provide little 
support for the hypothesis that this source is a ULX. In 
a study of 154 ULXs, the average luminosity is much 

lower (£ . 5 -8keV = 2 - 2 ± 4 - 5 x 1q39 er & s_1 )- The 
most luminous UL X ever detected, ESO 243-49 HLX-1 
(jFarrell et al.ll2009f ). has a 0.2-10 keV luminosity of 10 42 
erg s _1 , however, its hard X-ray (2-10 keV) luminosity 
is only 4x 10 40 erg s _1 because it has a soft photon index 
of 3.4. Therefore, if Mrk 739W is a ULX, it is the most 
luminous ULX in the hard X-rays (L2-10 kev) by over an 
order of magnitude. The average ULX power law index 
is also r= 1.97=1=0. 11, which is significantly softer than 
Mrk 739W (r=1.0 ± 0.2). Finally, no ULXs have shown 
such high amplitude variability over the short timescale 
of hours with this level of variability on ly seen on time 
scales of days to weeks (Gladstone 20l3). 

The measured SFRs in Mrk 739W provide an addi- 
tional constraint as to whether the hard X-ray emission 
could be from star formation. A SFR greater than 200 
Mq yr _1 would be required to generate the observed 
hard X-ray luminosity bas ed on the relationsh ip between 
SFR and X-ray emission (Ra nalli et al .1120031 ). The pre- 
dicted SFR in Mrk 739W is 0.3, 0.6, and 3.1 M© yr" 1 
from the Ha, UV luminosity, and 1.4 GHz emission. The 
predicted SFR from FIR emission of the combined sys- 
tem is 6.9 Mq yr _1 . All of these rates are significantly 
lower than the 200 M yr _1 needed to generate the ob- 
served hard X-ray luminosity. In addition, it is likely 
that much of this star formation would be extended and 
resolved in Chandra. 

It is interesting that Mrk 739W has not been detected 
as an AGN using o ptical emission line spectroscopy. 
iNoguchi et al.l ([2010) found that optical emission line 
studies are biased against 'buried AGN' that have a small 
scattering fraction or a small amount of narrow line re- 
gion gas. AGN with a low ratio of [O III] to hard X-ray 
luminosity (L[om]/^2-io keV <0.1) tend to be 'buried 
AGN'. The Z/[oni]/-^2-iokeV=0.008, consistent with a 
'buried AGN' and the lowest ratios found in their study. 
This finding is also consistent with a recent study that 
found that merging AGN selected in the ultra hard X- 
rays tend to have low L[ ni]/£i4-i95 keV ratios and be 
preferentially mis classified using optical line diagnostics 
(jKoss et alJEOlOh . 

For AGN with low luminosity in the [O III] line, 
nebular emission from star formation can overwhelm 
the AGN signature in o ptical emission line diagnostics. 
iSchawinski et al.l (|2010l ) found that for I/[ O ni]=10 40 erg 
s _1 , nearly 54% of star forming galaxies with AGN will 
be classified as star forming or composites. The small 
value of £[oiii]=7.5 x 10 39 erg s _1 in Mrk 739W suggests 
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Figure 4. Upper panel: Chandra spectra of Mrk 739E (left) Mrk 739W (right). Middle panel: Chandra count rate in 1.5 ks bins 
for Mrk 739E (left) Mrk 739W (right). There is a factor of 2.5 variability in flux over the 4-hour observation of Mrk 739W. Lower panel: 
Swift BAT spectra (left) and average count rate in 3 month bins (right) for Mrk 739. 



that star formation is overwhelming the AGN photoion- 
ization signature. 

3.5. CO Properties and Molecular Gas Mass 

CO velocity profiles can provide information on the 
dynamics of the molecular gas. The 12 CO 3-2 and 2-1 
spectra in Mrk 739 have almost identical shapes. Each 
spectrum has a narrow profile with FWHM=94±8 and 
98±6 km s _1 for single-Gaussian fits to the CO (2-1) and 
CO (3-2) profiles, respectively. These profiles are signif- 
icantly narrower than the CO (2-1) emission from NGC 



6240 (Figure 5) and imply a nearly face-on orientation 
to any disk-like structure in this system. 

CO measurements also provide estimates of the 
amount of molecular gas. In Mrk 739, Jco=109±33 and 
169±51 Jy km s -1 fo r the 2-1 and 3-2 lines. Following 
ISolomon et all (fl992h . L^ o =10.9xl0 8 and 7.5xl0 8 K 
km s _1 pc -2 for the 2-1 and 3-2 lines. Adopting a=1.5-4 
Mq (K km pc -2 ) -1 for the conversion from CO luminos- 
ity to molecular hydrogen, we find log; (Ho ) =9.2-9.6, 
similar to the Milky Way (jSanders et al.lll984f ). 
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3.6. Relative Velocity of the Binary AGN 

Measurements of radial velocities provide important 
insights about the dynamics of the merger. We use the 
Na I AA 5890, 5896 (Na D) absorption lines from stars 
and cold gas sin ce narrow emission l ines in AGN often 
have blueshifts ([Bertram et al.l [20071 ). There is an off- 
set of ~40 km s _1 between the two bulge components 
(8995±15 km s" 1 for Mrk 739E and 8953±15 km s" 1 
for Mrk 739W). The CO data also show evidence of two 
components with the peak brightness temperatures sim- 
ilar to the radial velocities in the Na D absorption lines 
(Figure 5). When fit with gaussians, the peaks consistent 
with the Na D radial velocities (8921±22 and 8980±16 
km s- 1 in CO 2-1 and 8956±12 and 8993±22 km s" 1 ). 
High resolution (1") inter ferometric CO imaging of this 
system would provide evidence to confirm this picture. 

There is also evidence of outflows in the narrow line 
region of Mrk 739E. There is a 192±22 km s" 1 blueshift 
in the [O III] line and a 153±25 km s _1 in the lower 
ionization [O I] A 6300 line compared to the Na D ab- 
sorption. This blueshift is consistent with other nearby 
QSOs which have an average [O III] blueshift of -174 
km s _1 ([Bertram et al.ll2007D . In Mrk 739W, there is no 
evidence of outflows in the narrow line region. 

3.7. Bolometric Luminosity and Eddington Ratios 

Using Chandr a and UV photometry and following 
IVasudevan et al.l (|2009l ). the bolometric luminosity is 
l.OxlO 45 erg s" 1 in Mrk 0739E (Figure 5). The ex- 
tinction corrected 2500 A luminosity is logL ° = 

J to 2500A 

43.7 =b 0.3. Using H/3 an d continuum emission 
([Vestergaard fc Peterson! l2006f ). the black hole mass 
is log Mbh=7.04±0.4, giving an Eddington ratio of 
AEdd=0.71, t he highest amongs all the Swift BAT se- 
lected AGN (Vasudevan et al. 2010). Our estimates are 
consistent with lHo et al.l (j2008[ ) who find AEdd=0.78 us- 
ing only optical spectra and the same method to deter- 
mine black hole mass. Uncertainties in intrinsic dust red- 
dening, as well as the inclination angle and spectral hard- 
ening parameter in the accretion disk model can lower 
the Eddington ratio at most 58% to AEdd=0.30. 

In Mrk 739W, the bolometric luminosity is 2xl0 43 erg 
s _1 using only the hard X -ray data with a bolom etric 
correction factor of 22 from IVasudevan et al.l ([2009). 

4. DISCUSSION 

We discovered a binary AGN in the galaxy Mrk 739 
based on Chandra imaging showing two unresolved 
(FWHM^300 pc) luminous hard X-ray sources with a 
projected separation of 3.4 kpc (5.8^0.1") . We find that 
a high level of star formation combined with a very low 
^[oni]/^2-iokeV ratio cause the AGN to be missed in 
optical spectroscopy. In the radio, there is resolved emis- 
sion with a spectral index consistent with star formation. 
The CO observations of the (3-2) and (2-1) lines indicate 
large amounts of molecular gas in the system. This gas 
could be driven towards the black holes during the violent 
galaxy collision and be key to fueling the binary AGN. In 
Mrk 739E, there is a high Eddington ratio (A E dd=0.71) 
and small black hole (logMeH = 7.05±0.3) consistent 
with an AGN accreting at a high accretion rate. Other 
than NGC 6240, this stands as the clearest and nearest 
case of a binary AGN discovered to date. 



Mrk 739 is an important example of how critical high 
resolution (< 1"), hard X-ray (>2 keV) imaging is in 
finding binary AGN (< 5 kpc). Observations with 
Chandra provide one of the most effective tools since 
obscuration and/or contamination from merger induced 
star formation can hide the AGN at other wavelengths. 
Mrk 739W showed no evidence for hosting an AGN un- 
til the Chandra observation despite a host of previous 
observations including UV and optical spectroscopy, and 
radio data from the VLA. While mega surveys such as the 
SDSS, 2dF, and 6dF are finding valuable information on 
hundreds of thousands of AGN using optical emission line 
diagnostics, this technique can be biased against finding 
AGN in objects that have high levels of star formation 
or obscuration like mergers (see also iKoss "eTall 1201(1 
IGoulding fc Alexanderil200l iVeilleux et al.ll2009h . 

The three nearest binary AGN (NGC 6240, Mrk 739, 
Mrk 463) discovered to date with Chandra share many 
properties which may hold clues as to why they form. 
A surprising result is that all three are luminous in the 
ultra hard X-rays and detected in the Swift BAT all sky 
survey (log £14-195 keV > 43.4 erg s" 1 ). This suggests 
that binary AGN happen more often in systems with 
bright X-ray AGN. More X-ray follow-up work needs to 
be done with less luminous merging systems to confirm 
this result. The large FIR (60 and 100 jam) luminosities 
(log L ,fi R =11.5, 11.1, and 10.6 for NGC 6240, Mrk 463, 
and Mrk 739) suggest that these systems may be linked 
to gas-rich progenitor galaxies consistent with theoretical 
models (Yu et al. 2011, submitted). 

Facilities: CXO, Swift, Sloan, IRAS, Gemini:GMOS, 
JCMT, VLA Facilities: Swift, Sloan, KPNO:2.1m, 
IRAS, AKARI 
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Figure 5. Upper left: CO observation of Mrk 739 with an observation of NGC 6240 overlai d for compa rison. The brightness temperature 
of NGC 6240 is reduced by 0.63 to account for the increased distance of Mrk 739 (Solomon et al. 1992). Upper right: CO observation 
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extinction-corrected SED model is shown as a solid line and the dashed line indicates the fit to the observed data. 
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